Thermal equilibrium and statistical thermometers in special relativity 
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There is an intense debate in the recent literature about the correct generalization of Maxwell's 
velocity distribution in special relativity. The most frequently discussed candidate distributions 
include the Juttner function as well as modifications thereof. Here, we report results from fully 
relativistic one-dimensional (ID) molecular dynamics (MD) simulations that resolve the ambiguity. 
The numerical evidence unequivocally favors the Jiittner distribution. Moreover, our simulations 
illustrate that the concept of 'thermal equilibrium' extends naturally to special relativity only if a 
many-particle system is spatially confined. They make evident that 'temperature' can be statistically 
defined and measured in an observer frame independent way. 
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At the beginning of the last century it was commonly 
accepted that the one-particle velocity distribution of a 
dilute gas in equilibrium is described by the Maxwellian 
probability density function (PDF) 



fM{v;m,P) = [ -^j cxp 



(1) 



[m is the rest mass of a gas particle, w e K'' the veloc- 
ity, T — (fce/?)"^ the temperature, fee the Boltzmann 
constant, d the space dimension; throughout, we adopt 
natural units such that the speed of light c — l\. When 
Einstein P, 0| had formulated the theory of special rel- 
ativity (SR) in 1905, Planck and others noted immedi- 
ately that /m is in conflict with the fundamental rela- 
tivistic postulate that velocities cannot exceed the light 
speed c. A first solution to this problem was put forward 
by Jiittner 0] . Starting from a maximum entropy princi- 
ple, he proposed the following relativistic generalization 
of Maxwell's PDF: 

/j(i;;m,/3j) = — 7(i;)2+'* exp[-/3jm7(t;)], \v\ < 1 

(2) 

[Z] = Z:]{m, f3j,d) is the normalization constant, E — 
m'y{v) = {m? + p^Y^'^ the relativistic particle energy, 
p = ■mv"f{v) the momentum with Lorentz factor "f{v) — 
(1 — v^)^^/"^]. Jiittner's distribution ([2]) became widely 
accepted among theorists during the first three quarters 
of the 20th century 0, H, @, H H - although a rigor- 
ous microscopic derivation is lacking due to the diffi- 
culty of formulating a relativistically consistent Hamil- 
ton mechanics of interacting particles [13, U, 12, T^ . 
Doubts about the Jiittner function /.j began to arise in 
the 1980s, when Horwitz et al. 13, HI proposed a 'man- 
ifestly covariant' relativistic Boltzmann equation, whose 
stationary solution differs from Eq. ([2|) and, in partic- 
ular, predicts a different mean energy-temperature rela- 
tion in the ultrarelativistic limit T ^ 00 [ld\. Since then. 



partially conflicting results and proposals from other au- 
thors [13, [18, 19, 20, 21] have led to an increasing confu- 
sion as to which distribution actually represents the cor- 
rect generalization of the Maxwellian ((T|). For example, a 
recently discussed alternative to Eq. ([2]) is the 'modified' 
Jiittner function [isl, 



to'* 7(1') 



2+d 



exp[-^Mj'7i7(t>)]. (3) 



The distribution ([3]) can be obtained e.g. by combining 
a maximum relative entropy principle and Lorentz sym- 
metry [20|. Compared with /j at the same parameter 
values /3j — /3mj ^ the modified PDF /mj exhibits 
a significantly lower particle population in the high en- 
ergy tail because of the additional l/£'-prefactor. 

Identifying the correct relativistic equilibrium velocity 
distribution is essential for the proper interpretation of 
present and future experiments in high energy and astro- 
physics 23, 23 , 24 , 2^ . Examples include the application 
of relativistic Langevin equations [26 



271 . |28l| to heavy ion 



collision experiments [22, 25], thermalization processes 
in ultra-relativistic plasma bearns |24l] , or the relativistic 
Sunyaev-Zel'dovich (SZ) effect [23], describing the dis- 
tortion of the cosmic microwave background (CMB) ra- 
diation spectrum due to the interaction of CMB photons 
with hot electrons in clusters of galaxies [2^ 113, [sH . The 
predicted strength of these spectral distortions and the 
cosmological parameters inferred from the SZ effect de- 
pend sensitively on the assumed electron velocity distri- 
bution [23]. 

Relativistic MD simulations.- To resolve the uncer- 
tainty about the relativistic equilibrium velocity PDF, we 
performed fully relativistic one-dimensional (ID) molec- 
ular dynamics (MD) simulations. The restriction to the 
ID case is inevitable if one wants to treat localized parti- 
cle interactions in a relativistically consistent manner, 
cf. remarks below. In our computer experiments we 
simulated the dynamics of classical, impenetrable point- 



2 



particles with elastic point-like binary collisions, employ- 
ing an algorithm similar to those of Alder and Wain- 
wright [s^l and Masoliver and Marro [s^. The basic 
time step of the algorithm involves three partial tasks: 
(«) determine the next collision event {xc,tc); (ii) evolve 
the system up to time tc', (in) calculate the momenta 
after the collision. The third task is solved as follows: 
If two particles A and B meet at the space-time point 
(xc, tc), then they exchange momentum according to the 
relativistic energy momentum conservation laws, 



Pa+Pb =PA+ Pb 
E{mA,PA) + E{mB,PB) = E{mA,PA) + E{mB,PB) 



(4) 



Here p = •mv^{v) is the relativistic momentum, and 
E{m,p) = (m^ + p'^y/^ the energy; hat-symbols de- 
note quantities after the collision. Given the momenta 
{pa,Pb) before the collision, the conservation laws (|4]) de- 
termine the momenta {pa,Pb) after the collision by [l9l | 

PA 'y{vof[2voE{mA,PA) - (1 + Vq)pa], 
Pb = ^{vQy[2vQE{mB,PB) - (1 + Vo)pb], 

where vq = {pA + PB)/[E{mA,PA) + E{mB,PB)] is the 
collision-invariant, relativistic center-of-mass velocity of 
the two particles. By assuming strictly localized, point- 
like pair interactions, one may avoid the introduction 
of fields which are required when considering relativistic 
particle interactions-at-a-distance (the interested reader 
may wish to consult the original papers of Wheeler 
and Feynman Currie et al. [13], and Van Dam and 



Wigner pjJ,|3J], who discuss in detail the difficulties asso- 
ciated with classical particle-particle interactions in SR) . 
However, considering point-like localized interactions is 
expedient in the ID case only; in higher space dimen- 
sion the collision probability would become zero, thus 
preventing the system from equilibration. Moreover, if 
two colliding particles carry the same rest masses then 
elastic ID collisions merely interchange their velocities; 
hence, elastic binary collisions are not able to drive a 
ID one-component gas to equilibrium. In our simula- 
tions we considered a two- component mixture, consisting 
of A^i light particles having equal masses mi, and N2 
particles with equal masses 7712 > mi . The motion of the 
N = Ni -\- N2 particles was restricted to the ID inter- 
val [0, L], assumed to be stationary in the lab frame S. 
The results presented below refer to elastic reflections at 
the boundaries; however, we found that periodic bound- 
ary conditions yield identical outcomes if the total initial 
momentum was chosen to be zero in E. Generally, our 
simulations mimic a relativistic microcanonical ensem- 
ble, since the total initial energy i?tot in E is conserved 
in the microscopic collision processes. The above conven- 
tions define the simplest interacting model system that 
{i) complies with all principles of SR, (ii) does not require 
the introduction of interaction fields, (Hi) can be simu- 
lated without further approximation, and (iv) exhibits a 



universal stationary equilibrium state. Hence, this model 
systems provides an optimal test case for probing the pre- 
dictions of different relativistic kinetic theories by means 
of numerical experiments 0, 0, 0, H, 0, Moreover, 
as we shall see below, it helps to clarify longstanding con- 
troversial questions regarding the definition and meaning 
of 'temperature' (i.e. thermometers) and 'thermal equi- 
librium' in SR. 

Numerical results.- In order to identify the station- 
ary one-particle velocity distributions for the light and 
heavy particles, respectively, we waited until the ID two- 
component gas had approached the equilibrium state 
(typically, after 10^ collisions per particle). Then the 
particle velocities were measured E-simultaneously, i.e., 
at equal times with respect to the resting lab frame E. 
To increase the sample size we repeated this procedure 
several times during a simulation run and collected the 
data into a single histogram. An example is shown in 
Fig. [U based on a simulation with N = 10000 parti- 
cles (A^i = N2 — 5000, m2 = 2mi). Each particle had 
been given a random initial position Xi(0) G [0, i] and 
a random initial velocity Vi{0) = ±0.8, corresponding to 
a mean energy per particle e = 2.5mi. As evident from 
Fig. [H for both particle species the numerically obtained 
one-particle PDFs (o) are in very good agreement with 
the standard Jiittner function /j (solid line), and differ 
significantly from the modified distribution /mj (dashed 
lines). The same result was found for A^i ^ A^2- 
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FIG. 1: Equilibrium PDFs in the lab frame E: Numerically 
obtained one-particle velocity PDFs (o) based on a simulation 
with A''i = 5000 light particles of mass mi and N-^ = 5000 
heavy particles with mass 7712 = 2mi. The mean energy per 
particle in E is e = Etot/{Ni + N2) = 2.5mic^ The solid 
curves correspond to Jiittner functions ^ with same pa- 
rameter /3j — 0.702 (mic^)~^, but different particle masses, 
respectively. Dashed lines show the corresponding modified 
distribution Q with /3mj = 0.402 (mic^)"^. As the distri- 
butions are symmetric with respect to the origin, only the 
positive velocity axis is shown. The simulation data is con- 
sistent with the standard Jiittner distribution and thus 
provides evidence against the modified distribution (|3j. 
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The distribution parameters Pj/mj were determined 
from the initial energy by means of the foUowing con- 
sideration: If the particle numbers A^i and N2 are 
sufhciently large (thermodynamic limit), then the one- 
particle PDFs in the lab frame S are expected to converge 
to either /j from Eq. ^ or /m j from Eq. ^ . Generally, 
the mean relativistic energy value /i of a one-particle PDF 
f{v; m, (3) is given by 



d'^v f{v; m, f3) mj{v). (6) 



{kl<i} 



Assuming (i) that an equilibrium state exists where both 
species can be described by the same value (3, and (ii) 
that for a gas in equilibrium the mean energy per particle 
is the same for particles of the same species, the total 
energy can be expressed as 



Etot = Ni fi{mi,l3)+N2 A*(™2,/3)- 



(7) 



In our case, the energy mean values of the two ID can- 
didate PDFs /j and /mj read explicitly 



^mj(to,/3mj) = m 



2ifi(/3jm) 
Ki{l3Mjm) 
^o(/3MJ"^)' 



(8) 



with Kn denoting the modified Bessel function of the 
second kind [35| . For each simulation run the parame- 
ter tuple {Etot, Ni, N2,rni, 1712) is known. Hence, upon 
inserting them into Eqs. ([7]) and ([8]), these parameters 
uniquely determine the parameter value Pj/mj that is 
consistent with the chosen velocity PDF /j /mj ■ 

Temperature and equilibrium.- Most remarkably, in 
spite of the different particle masses the two numerically 
obtained velocity PDFs in Fig. [1] are very well matched 
by Jiittner functions ^ with the same parameter 
According to our simulations, this holds true with high 
accuracy for a wide range of initial conditions and mass 
ratios. Hence, the Jiittner function does not only pro- 
vide the best 'fit' to the numerical data, it also yields a 
well-defined concept of 'temperature' in SR: Intuitively, 
the temperature T is thought to be an intensive quan- 
tity that equilibrates to a common value if two or more 
systems are brought into contact with each other (i.e., 
may exchange different forms of energy). In our case, it 
is natural to consider the particle species as two different 
subsystems that may exchange energy via elastic collision 
processes. After a certain relaxation time, the combined 
system approaches a 'thermodynamic equilibrium state', 
where each subsystem is described by the same asymp- 
totic, two-parametric velocity PDF /j(z;; m^, /?j), differ- 
ing only via the rest masses rrii. The commonly shared 
distribution parameter /?j may thus be used to define 
a relativistic equilibrium temperature T :— (fce/^j)^^- 
However, for this concept to be meaningful, a restric- 
tion of the accessible spatial volume is required - be it 



by means of periodic boundary conditions, or by impos- 
ing reflecting walls. Otherwise, it cannot be expected 
that a many-particle system approaches a universal sta- 
tionary state which is independent of the specific initial 
conditions. This observation has an important impli- 
cation: Any (relativistic or non-relativistic) Boltzmann- 
type equation [1, H, [11 [H, [H, [H, 113 that gives rise to 
a universal stationary velocity PDF implicitly assumes 
the presence of a spatial confinement, thus singling out a 
preferred frame of reference. 

Moving observers and statistical thermometers. - 
From our simulations we may further determine the equi- 
librium velocity distributions as seen from another frame 
S' moving with velocity u relative to the lab frame E. 
Figure [2] depicts the results for u — 0.25 and same simu- 
lation parameters as in Fig. [TJ In contrast to Fig. [2 the 
numerical data points in Fig. [2] were obtained by mea- 
suring velocities E'-simultaneously. The solid curves in 
Fig. [2] correspond to the PDF 

f,{v';m,(3j,u) = exp[-f3Mn)mj{v') (1 + W)] 

Zj 7(m) 

(9) 

[v' is the particle velocity in the moving frame E']. The 
PDF Q reduces to the Jiittner function ([2]) for u = 0; 
/j is obtained by using the fact that the one-particle 
phase space PDF, reading 

$j(x,p) = (ZjL)-i exp[-ATS(m,p)] e(x) Q{L - x) 

in E, is a Lorentz scalar [sl, [s^; 8 is the Heaviside unit- 
step function. Due to the excellent agreement between 
the numerical simulations and Eq. ([9]), we may state more 
precisely: Two relativistic gas components are in 'ther- 
modynamic equilibrium' for any observer if their one- 
particle velocity PDFs are given by generalized Jiittner 
functions ^ with same parameters /3j and u. Only in 
this case the net energy transfer between the different gas 
components in the container vanishes. Last but not least, 
the above results shed light on a longstanding, highly de- 
bated question 43j.41, 42, 4^ 44, 4^ originally posed by 
P. T. LandsberggJ"^ 

Does a moving body appear cool? - Evidently, the 
answer depends on the thermometers employed by differ- 
ent observers. Adopting, for the reasons discussed above, 
T :— (fcB/3j)~^ as a reasonable temperature definition, a 
moving observer with rest frame E' can measure T by ex- 
ploiting the Lorentz invariant equipartition theorem [4^ 



(10) 



where u — —{v')', and averages ( • )' are taken E'- 
simultaneously. We verified the validity of Eq. (|10p ex- 
plicitly by using simulation data obtained for different 
values of u. Hence, Eq. (jTU]) defines a Lorentz invari- 
ant gas thermometer on a purely microscopic basis. Put 
differently, this intrinsic statistical thermometer deter- 
mines the proper temperature of the gas by making use 
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FIG. 2: Equilibrium PDFs in a moving frame E': Velocity 
PDFs as measured by an observer who moves with veloc- 
ity u = 0.25c relative to the lab frame E. Parameter val- 
ues are the same as in Fig. [T] The solid lines correspond to 
Jiittner functions /j from Eq. (|9]) with the same parameter 
/3j = 0.702 (mic^)-! as in Fig. [T] and different masses mi and 
m2, respectively. 



of simultaneously measured particle velocities only; thus, 
moving bodies appear neither hotter nor colder. Analo- 
gous considerations apply to the 2D/3D case. 

Summary.- Fully relativistic MD simulations favor 
the Jiittner distribution ^ as the correct relativistic 
one-particle equilibrium velocity distribution. The re- 
sults are conclusive for the ID case, and provide evi- 
dence against theories [13, E El, [H H Q that pre- 
dict other distributions. Further, our simulations corrob- 

that the tempera- 



orate Landsberg's hypothesis [42, |46 
ture of classical gaseous systems can be defined and mea- 
sured in a Lorentz invariant way. The extension of the 
MD approach to higher space dimensions is nontrivial, 
due to the fundamental difficulty of treating 2D and 3D 
two-body collisions in a relativistically consistent man- 
ner [lOi [ul [ij, 111]. As a first step, it should be 
carefully analyzed if and how specific semi-relativistic in- 
teraction models affect the 2D/3D equilibrium velocity 
distribution. 
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